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| By Donald F. Hays 
Research Laboratories 


A Composite Infinite Slider Bearing 





General Motors Corporation 


An analysis is made of an infinite slider bearine which is com- 
posed of three intersecting surfaces. Characte.istics of such a 
composite slider are then compared with the corresponding charac - 
teristics of an infinite plane slider with the same film ratio. 
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Film height at surface intercepts 
Friction force - F/L 
Integration constants 
Slider length - L 
Pressure - F/L* 
Relative velocity - L/T 
Load Capacity - F/L 
Center of pressure - L 
Coefficient of friction 
Film Thickness - L 
Exit film thickness - L 
Inlet film thickness - L 
Film ratio = h,/h, 
Subscript - 1, 2, 3 range 
Flow Rate - L’/T 
Velocity - L/T 
Coordinate axes 
Amplitude ratio = 5/h, 
Transformed coordinates 
Viscosity - (FT)/L* 
Angle of inclination 
Waviness amplitude - L 


Note: Circumflexed symbols are dimensional ( ) 














DONALD F. HAYS 
INTRODUCTION 


In the manufacture of bearings, care must be exercised in the 
preparation of the bearing surface. Two factors, surface roughness 
and surface waviness, must be taken into consideration. Surface 
roughness may be considered as a high frequency irregularity of the 
surface in contrast to waviness which is a low frequency phenome- 
non. Since waviness is sometimes found in ground bearing sur- 
faces, it is of interest to investigate the effect of this irregularity 
on bearing performance. This is especially important when the 
amplitude of the waviness approaches an order of magnitude which h 
is comparable to the minimum film thickness. The simplified model 
to be analyzed is an infinite composite slider bearing which is com- - 
posed of three surfaces. These are orientated such that they form a 
linear wave form about a plane surface which is taken as the refer- 
ence slider. For various amplitude ratios based on the outlet film 
height, the operational characteristics of the bearing are computed 
and are compared to the corresponding characteristics of the infi- 
nite plane slider. In all comparisons, the ratio of inlet film thick- 
ness to outlet film thickness is constant. The computed character - E 
istics include the pressure distribution, load capacity, frictional a 
force, coefficient of friction, center of pressure, and flow rate. The 
assumptions are that the fluid is incompressible, the viscosity is 
constant, and that the Reynolds’ Number is small. 





ANALYSIS 


The infinite composite slider is shown in Fig. 1. It is composed 
of three sections of length L/3, the sections being inclined at the 
angles 8,. The entrance film thickness is h, and the outlet film t 
thickness is ho. A separate coordinate system is used for each sec- 
tion and each section is solved with the appropriate boundary con- 
ditions sc as to insure a continuous pressure function throughout the 
total length. The basic equation which must be satisfied for each 
section is the one-dimensional Reynolds’ Equation. 


d 

ax, (Ys ax) * “one a a) 
The subscripts “n” have a range 1, 2, 3, and indicate the section be- 
ing considered as measured from the trailing edge. The film thick- 


ness is given by the following form, 


hy Bn Xn (2) 
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Figure 1. Diagram of Slider Bearings. 


Equation (1) may be integrated twice to yield the pressure equations 
as shown. 
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the nondimensional pressure may be written as 
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These Equations (5) contain six arbitrary constants of integration 
which are evaluated such that the following six boundary conditions 
are satisfied: 


6 =0 P,=0 (6) 


(7) 
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-1 6 =0 P =P, (10) 
1 @& =0 P, =P, (11) 


After considerable algebraic manipulation, the constants of integra- 
tion take the following form: 


Ki. * 2M (12) 
K,,° -G . 1) (13) 
K 2 > (14) 
K,, = 3 (™ - 4) (16) 
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where M and N are given by the following expressions: 
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The expressions for A, and B,, are determined by the geometry of 
the slider as shown in Fig. 1. By defining a nondimensional ampli- 
tude of waviness as 


my 
A= he (20) 
and the film ratio as 
K = = (21 
ho 


then the expressions for A, and B, become 





A, =1 (22) 
1 

A, = 3 (2 + K + 34) (23) 
1 

A, = 3 + 2K - 3A) (24) 

B, =K-1+ 34 (25) 

B,=K-1 - 64 (26) 

B, = K - 1 + 3A = B, (27) 





The load capacity of the composite slider bearing is found by 
summing the load capacities of each section. Thus 
2 ) 
w= 2 f Pax, (28) 


n=1 


) 
} 
LOAD CAPACITY 
} 
| 
| which when written in nondimensional form becomes 
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DONALD F.HAYS 
FRICTIONAL FORCE 


The frictional drag experienced by the bearing as it moves with 
a velocity *U” is given by 


I 





Ss 5 du 
j n 
F = - os lu J, ry dX, (30) 
where 
- 1 dP, r ( y 31 
Un = 3a ax, 9 - nly - U “_ (31) 


In nondimensional form the friction force per unit width is 











38, Kin Pn 
uP 4in(1 + 32 
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CENTER OF PRESSURE 





The center of pressure of the composite slider as measured 
from the exit section is found from the following expression: 





3 
5 (n-1)L 
an J, P15, sa 3 ax 
X = (33) 
\* ¥ 
uf Pz dX, 


The denominator of Equation (33) is the load capacity as given by 
Equation (29) and the numerator of Equation (33) is given in nondi- 
mensional form as follows: 





( ) ho 6 dD 1 ({n-1} 3 2) 
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~N 


The nondimensional center of pressure as measured from the trail- 
ing edge is 


X _ Equation (34) (35) 
L Equation (29) 


COEFFICIENT OF FRICTION 





The coefficient of friction defined as the ratio of frictional re- 
sistance to normal load is given by 


L ,? _ Equation (32) 
ne Equation (29) (36) 


FLOW RATE 


The flow rate per unit width is determined from 








- -hy dP, Uh, 37 
1” i2y aX, a (37) 
xX) to) 
which reduces to the following nondimensional form 
@ «By 
S ae 
Uh, 5 (38) 
RESULTS 


The characteristics of the composite slider bearing have been 
computed for film ratios up to 5.0 and for various amplitude ra- 
tios A. Since these amplitude ratios may be either positive or nega- 
tive, two curve sets are given for each characteristic of the slider. 
It should be noted that for negative values of A there is a limiting 
value for k such that when 


K< 1-34 (39) 


there is an internal film thickness that is less than the outlet film 
thickness. In all curves, the ratio of the composite slider charac- 
teristics to the plane slider characteristics is made on the basis of 
the film ratio which is computed from the entrance and exit film 
thickness. For positive values of A the limiting case is for 
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K Sat] (40) 


which also yields an internal film thickness that is less than the exit 
height. 

In Figs. 2 and 3 are shown pressure curves for a slider with a 
film ratio of 3.0 and having various amplitude ratios A. The opera- 
tional characteristics of a composite slider as compared with an 
ordinary plane slider are shown in Figs. 4 through 13. In Figs. 4 
through 13, “Region I” indicates a minimum film thickness occurring 
at the exit section and “Region II” indicates an internal minimum 
film thickness. A symbol such as W(o) refers to the characteristic 
of the ordinary plane slider, i.e., A is zero. 


CONCLUSIONS 





The overall effect of waviness in the surface of an infinite plane 
slider is seen to be deleterious to the operation of the bearing. This 
effect becomes more pronounced as the amplitude ratio is increased 
for a value of the film ratio. Thus waviness present in the surface 
might not be detectable in its influence on the bearing operation at 
light loads, but could possibly contribute greatly to the reduction in 
load capacity at higher loads when the minimum film thickness de- 
creases. The general effect of waviness is to aggravate the situa- 
tion; higher loads decreasing the film thickness, a larger amplitude 
ratio then decreasing the load capacity of the unit which in turn 
would decrease the minimum film thickness. This cycle would con- 
tinue until equilibrium was established. 

For extremely thin films, the viscosity of the oil will continue to 
decrease due to the effect of temperature as the oil passes through 
the bearing. Since this analysis is based upon a constant average 
viscosity, some discretion must be exercised in the interpretation 
of the results. 
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Geometries for Positive Amplitude Ratios. 
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A Two-Dimensional Theory of Stress 
Concentration Around Voids 














By Leonard G. Johnson 
Research Laboratories 
General Motors Corporation 


INTRODUCTION 





Present day emphasis on the various aspects of product relia- 
bility has led, in particular, to a serious consideration of stress 
raisers (holes, inclusions, etc.) in engineering materials. These 
stress raisers can very well be one of the chief causes of early 
failures in service. Therefore it behooves the materials engineer 
to become familiar with the magnitude of stress concentration pro- 
duced by holes and inclusions of various shapes and properties. The 
present investigation will be confined to the problem of determining 
the Cartesian components of elastic stress around a hole in a plate. 


STATEMENT OF THE PROBLEM 





We assume a plate is subjected to a cuustant two-dimensional 
stress system whose Cartesian components are 


(1) S, = Normal stress in x-direction (horizontal) (1) 
(2) S, = Normal stress in y-direction (vertical) (2) 
(3) T,, = Shear stress component (3) 


We represent this by the type of sketch shown in Figure 1. At great 
distances from the hole the components of stress are (S, , S,, T,,,,). 
However, in the neighborhood of the hole the stress components will 
be considerably different from (S,, S,, T,,,). In general, let us as- 
sume the components of stress at any point (x,y) are (7,, Gy, Tay)- 
Then if the origin of our coordinate system is assumed to be at the 
centroid of the hole, we can state that 


lim 9, = 8, (4) 
r= « 

. 2 2 

lim 9y =Sy (r = Vx’ +y’) (5) 
re © 

lim ty = Txy (6) 
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Let Figure 2 represent an enlarged view of the hole. Suppose the 
equation of the boundary of the hole is y = Rix). We are assuming 
a no-load boundary such that for each (x,y) on y = B(x) 


Radial stress = 9. = 0 (7) 
Shear stress = T.g = 0 (8) 
Converting these into Cartesian components (Rp % ,p°%,p%,, ) on 


the boundary of the hole we obtain the following conditions along 
y = B(x): 


B% /B% = [B’(x)}" (9) 
BYx* BY, B T xy (10) 


STRESS FUNCTION NOTATION 





It is well known |Ref. 1] that the Cartesian components of elastic 
stress can be represented by the second partial derivatives of a 
function ® which satisfies boundary conditions and the biharmonic 
partial differential equation 


. 2 am. a = 
axt * 2 aeF ay? pyr = 9 (11) 
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Figure 2 


Then the stress components are 








y 

; r,t 
} ai dy* (12) 

} 820 
% °* “SS (Assuming no body forces) (13) 

| - a? 
o’xy ~ “Ox ay ates 


If we let (x,y) be the stress function for the plate without the hole 
| it follows that (12), (13), and (14) must represent S,, S,, and T,,, 
respectively. 
Now let us assume that y(x,y) = (x,y) + A(x,y) is the stress 
function for the plate with the hole. Then the stress components be- 
} come 


02a 
+ £5 (15) 
dy 
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a*y a° aA qv A 
- * ax* axe ox Sy * 73x (16) 
2y 26 2x =r 
"ay © « Ro ma : - : = Tky - amr (17) 
; ox ay oxay oxay ' axay 


This form is assumed in order to retain the Cartesian frame of 
reference throughout the investigation. Hence, on the boundary of 
the hole (along y = B(x)), the conditions (9) and (10) become 





ax: 2 
a= (B’( (18) 





a7A O7A\2 
> ~ @x ay (19) 





Although A is not readily discovered in any particular problem, this 
form of solution does indicate one possible avenue of approach. 


SPECIAL CASE OF A CIRCULAR HOLE 





For a circular hole of radius a > 0, the equation of the boundary 


r+; ce (20) 
For the locus (20), 


dy B’ (x) 2 


x 
ax -y (21) 


Using (21) in conditions (18) and (19), we find that we can take 











patel S = 7 x a 1 a? | 2.2 
AUGY - x* +y* - . (x*+y*)” , a’ In('+y ) 
a*y* 1 x" -y* '- 
4 Ss ——— 4 = o® aaa —  - 2 2 4 
tc 4? Gxtsy2)? 4 2 ine v)] 


[ 2a*xy a*xy | 
°F - : 
5s L xe? (x ry?) | (22) 





(S., S,, and T.. are constants.) 
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Substituting the partial derivatives of (22) into (15), (16), and (17) we 
obtain the stress components (a # 0, x* + y” > a*) 


[1 3X -Gxty*ty® |, xi -Oxiy’ty’ 31. * 
f = _ : - site a — 
J, Ss, + D) a (x? +y? )* ¥ (x? +y ») 22 (x? +y”) 
3 , x? .6x? y’ +y* : x? -6x* y? +y* 1 , x? -y* 
+ Sy | - r tors (x? +y?)* +a (x? +y?)° +,8 (x? +y?)? | 





| 12a* xy(x?-y?) 8a*xy(x? -y’) 4a*xy 
+ Txy (x2 +y?)* - (x? +y? Ny - (x?+y?)° 


x* -6x?y *+y* x*-6x*y *+ x?-y? 
, - 8, [ -2et tery ,, ee 1] 














2° (x+y?) (x?+y?) 2° (x*+y?)° 
r x* _6x 2y? +y4 x* _6x? y 2+y4 x2 _y2 
+Sy |! +2 gS | a 8 tar | 
y 2 (x? +y*) (x*+y *) 2 «(x?’+y7) 
-12a* xy(x?-y*) 8a* xy(x* -y*) 4a? xy 
* Ixy | (x 2+y?)* (x? +y?)” . (x? +y?)? (24) 
y= 8, | 6 o2y?-y?)  xye?-y*) gy 
lida (x?4y2)© ~ * (xPay2)?  * xtay?)® J 
4 Xy(x*-y 4) 2 Xy(x?-y*) — 
+S, - 6a “(x24y2)) + 4a (x?+y2)> -a (x?4y2)2 
x4 -6x*y* +y* x4 -6x*y’ +y* 
+ Tyy [1 - 3a* a + 2a? san | (25) 


On the boundary of the circular hole (where x* + y* = a*) these com- 
ponents simplify to 


=% @ (3-40) }+S, y [3 - 4¢ e% |- 8T xy OE (26) 
(a # 0) 
B% = 8G) (3-40) ]+8,Q (3-48) ]-8Ty@'® (27) 


(x* + y* =a’) 


x x,2 x “ 2 2 
B xy ~ - 8, C¥) [3 -4() | - 8, ) [3 -4) }+ 8Txy eC) oe (28) 
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FORMULAS FOR “EFFECTIVE STRESS” 
ON A VOID BOUNDARY 








For two-dimensional problems involving holes in very thin 
plates (plane stress) or very thick plates (plane strain), it is often 
of interest to know the value of the “effective stress” defined by 





y2 5, = 2 | -0,) + (0, -0,) + @,-9,)° 
+ 6 (72, + 72, + 72, ) | (29) 
Since for plane stress: %, = \&« = Ty, = 0, we obtain 
y2se, = (0x +O) + 3 (Ty - Fy Dy) (30) 
For plane strain, where %, = v(C, +%,) and %s = TW, = 0, (29) re- 
duces to (v = Poisson’s ratio) 
ett, = (L-v+v") x+y) +3 (Tey - 9% Oy) (31) 
On the void boundary BT, -B°B?,=0. (See (10).) 
Hence, 
For plane stress: p%.rs =|p% +p%! . (32) 
For plane strain: p%.55 =|p%, + p% Vi-vev. (33) 
CONCLUSION 
For any void, y = B(x), in a plate subjected to stresses S,, S,, 


T,, it is possible to select an auxiliary stress function A(x,y) which 
will satisfy (18) and (19). The second partial derivatives of A(x,y) 
when added to S,,S,, T,,, in accordance with (15), (16), and (17) 
will then give the Cartesian stress components anywhere around the 
void, and questions of stress concentration can be answered. 
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Stability Study of a Fluid Flow Problem 
; with Flow Reversal 


By Richard L. Justice 
} Research Laboratories 
General Motors Corporation 


In considering the problem of fluid flow an interesting mathe- 
matical problem arises when the flow reverses. In order to gaina 
} thorough understanding of this process a very simple two volume 
model will be studied. This model is designed to introduce into the 
mathematics only those problems of instability and not consider the 
; problem of heat transfer. Insight gained from this problem may 
well lead to applications in real problems involving mathematical 
instability. 





TABLE OF SYMBOLS 





- Cross Sectional Area 
Specific heat at constant volume 
Specific heat at constant pressure 
Diameter of tube 

- Friction factor of tube 
- Gravity constant 

- Ratio of specific heats 
System constants 
Mass 

Reynolds’ number 
Pressure 

Heat 

- Gas constant 

- Temperature 

- Gas velocity 

- Volume 

- Mass rate of flow 
Work 

6 - Time 

p ~- Density 

u - Viscosity 
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Figure 1 


Applying the first law of thermodynamics to the volumes as 
shown in Figure 1, the following equations are obtained: 
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Since the process is considered adiabatic, —_—* 0. Further, let 
dW, dv dW, 2 
W (flow rate of the medium), and a = pf, = and —— -=0 


(since V, is constant). At a given instant if the piston is moving to 
the right, equations (1) and (2) become 
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dP, 1 d v:) 
r Oo Vv, (. KR T+ Wour - © P: aa (3) 
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de 7 Vv KR Tin Win (4) 
3 
Observe that as the piston traverses its allowed path, several 
conditions with regard to the pressures may exist; namely, P, > P, 
or P, P, with either turbulent or laminar flow depending upon the 
magnitude of the inequality. For the moment let us assume that 
P, > P, and the inequality is sufficient to develop turbulent flow. 
Define the following dimensionless parameter. 
P, = Vi - P,/P, 
From this, then 
qP, 1 1 dP, =; <*.) 5 
dé = 2 P, (- P, dé Pp? dé (9) 
and 


Next observe that for turbulent flow 
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By substituting in (5) the following expression is obtained 
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d P, Bara (%, , pi) k (1 - P2) dV, 
da 20O 2V; V; "—— Cnmwa 
dM, 
Since M, + M,= My and —_ W;,, = Wout We May now sum- 
marize the three essential equations for P, > P, turbulent flow. 
dP, kRKVYT,P,P, KP dV, - 
Oo - V, ~“V, “a 
dP, ew ato é' 2) ka -P) dv, " 
oe "-—Tv, \v, * - Pi) -—sE-v we 
dM, K P, P, 
ao °° (8) 
VT, 


Next consider the case of turbulent flow for P, > P, . Define 
P, = V1-P,/P, 


and hence, 


The following development parallels the above. 
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The turbulent P, > P, equations are summarized here 
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These equations (6)-(8), and (9)-(11), will not be applicable in 
the laminar flow region, since they are inherently unstable at the 
flow reversal point. This is caused by the occurrence of P, in the 
ao *. (Note that the a equa - 
tion will not be unstable.) Hence a new definition of the dimension- 
} less P, is required. 

Consider iaminar flow with P, > P,. Define 





denominator of the second term in 


p.=1 +2 
2 ? P, 
Then 
dP, 1 dP, P,dP, 
—— &. S|) S ee (12) 
dé P, dé P, dé 
dP d P, ' , , 
where — and = are defined as in (3) and (4). According to such 
H 0 


references as McAdams! the following relationships hold, 


‘McAdams, W.H., “Heat Transmission” 3rd Ed., McGraw-Hill, 


New York, 1954. 
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Then flow rate is obtained as follows 
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The following three equations hold for P, > P, laminar flow. 
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The final set of equations needed then are the laminar equations 
for P, > P,. In this case the following relationships apply. 
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Figure 2 illustrates the application of these four sets of equa- 
tions and the points at which a decision must be made to switch sets. 

In order to arrive at a satisfactory numerical solution of these 
equations, it is necessary to maintain two different types of ac- 
curacy. Firstly, it is necessary to hold the error within reasonable 
engineering accuracy, such that results could be corroborated by 
tests. Secondly, it is necessary to maintain an accuracy necessary 
for stability of the numerical technique employed to solve the dif- 
ferential equations. This latter one is by far the more critical one 
and is such that the “roots” must be located quite accurately; roots 
are the points of the cycle where a transition from turbulent to 
laminar or laminar to turbulent flow occur. The principal obser- 
vation of this study is that if these roots are located with sufficient 
accuracy then it is not particularly difficult to maintain a decent 
engineering accuracy over the rest of the cycle. 

With regard to the numerical techniques and difficulties, the 
Runge -Kutta-Gill method was used with root iterations to find the 
above mentioned roots. (This problem was solved with the use of an 


EQUATIONS 
TURBULENT 6, 7, 8 
+W LAMINAR 13, 14, 15 
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Figure 3 


IBM 704). It should be noted that the location of the roots is not a 
function of the independent variable alone, but is also a function of 
the dependent variables. This can lead to the following problem. 
When RKG has produced a number which appears to have spanned a 
root, an effort is made to iterate and locate the root more accu- 
rately. By this more accurate technique the root is determined 
more accurately also and actually appears to move to a point later 
in time. This then causes the program to stop. Hence it appears to 
be necessary to start decreasing the step size well ahead of the oc- 
currence of any root. This is accomplished by testing at every time 
step the relative merit of halving or doubling the step size on the 
next step. This costs heavily in compute time, but allows the ma- 
chine to operate at nearly optimum step size at all times. 

One of the oddities of this problem is the ratio of time in lami- 
nar flow to time in turbulent flow. On a given standard run, see 
Figure 3, it was observed that through 86% of real time the flow is 
turbulent while only 14% is laminar. However, due to stability prob- 
lems, 57% of the computer time is spent in turbulent and 43% in 
laminar. The anamoly here is that in general the usefulness of such 
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a program as this is greatest in the region of greatest change in 
pressure and hence highest flow. Therefore, the contribution to 
useful results obtained from the laminar region is negligible even 
though it costs about half the computer time. 

This problem was formulated with the hope of gaining a better 
understanding of the mechanism of flow reversal. The primary pre- 
requisites of the program were to get some measure of the stability 
requirements and the length of time necessary to compute converged 
results. The problem serves primarily as an initial step towards 
handling n linked volumes in series with n - 2 constant volume 
tanks and two varying volumes. The final program of interest then 
will be to yield a complete pressure, temperature and mass distri- 
bution picture of such a system as a function of time. 
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Computers —Large or Small? 


By Robert K. Louden 

System Engineer 

International Business Machines Corporation 
Detroit, Michigan 


INTRODUCTION 





If you were running a department within a company and your de- 
partment required electronic data processing, would you attempt to 
share a large computer with other departments, or would you prefer 
to go it alone, with a smaller computer set up specifically for your 
own department? This paper is a qualitative discussion of some of 
the factors which affect the choice of “Computers — Large or 
Small?” based upon the author’s experiences with various compa- 
nies, some of which chose large computers and some of which chose 
small computers. The opinions below are the author’s and do not 
necessarily reflect the judgments of any particular computer user. 

Small computers are here defined as digital systems renting for 
less than $5,000.00 per month and large computers are defined as 
digital systems renting for more than $25,000.00 per month. It is 
assumed that the data processing work of the department in ques- 
tion can be handled on either system; if the small computer lacks 
either enough processing speed or memory capacity to do the work, 
then the decision must automatically go to the large machine. It is 
also assumed that both systems are typical of new equipment being 
installed in 1961; i.e., that the large computer is a transistorized 
magnetic core system with multiple tape channels and off-line input 
and output, and that the small computer is a transistorized magnetic 
core system with buffered input and output. 


Cost and Throughput 





In the area of direct operating cost, or answers per dollar of 
computer rental, the large computer generally comes out well ahead 
of its small competitor. This is due to the fact that both large and 
small computers must possess certain basic input/output and logi- 
cal capabilities, and a comparatively large proportion of the small 
computer cost must be set aside to meet these basic requirements. 
The additional cost of the large computer can then be used to pro- 
vide increased speed and memory capacity which will produce sig- 
nificantly higher answer per dollar ratios. 
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However, if the large computer will be located some distance 
away from the departments which will use it, the cost of transpor- 
tation of men and information to and from the computer, plus the 
cost of time delays caused by the transportation problems, may be 
sufficient to overcome the large computer’s direct cost advantage. 
The cost of transportation and time delays has been explored at 
some length by Donald C. Gause of the General Motors Truck & 
Coach Division of General Motors Corporation, who had developed 
comprehensive data on the breakeven cost of large versus small 
computers as a function of the large computer’s distance and trans- 
portation time, considering transportation methods ranging from 
men on foot to telephone lines. 

TELE-PROCESSING equipment which allows high speed long 
distance transmission of data provides a satisfactory means for the 
transfer of low volume high priority data, but at present most high 
volume Tele-Processing Systems cost about as much as a typical 
small computer system. 

The greater speed and capacity of the large computer will en- 
able it to handle big jobs which might develop in the future, while 
the small computer may not be able to do so, due to its more limited 
memory size. Therefore, the likelihood of such big jobs developing 
in the department should be ascertained, if possible, before a de- 
cision on computing equipment is made. 


Communications 





Computer communication refers to the task of instructing the 
computer how to solve problems (called programming) and then 
checking and testing the instructions until they produce correct an- 
swers (called debugging). A programming language or code is used 
to write the instructions in a form which the computer can translate 
into basic computer operations such as add, subtract, etc. 

Usually, the large computers will utilize languages which are 
very general and very powerful. The small computers will use lan- 
guages which are (by comparison) simple and ridigly defined. This 
usually means that experienced programmers working on large 
complex problems will have an easier time of it with the big ma- 
chine languages, but that novice programmers on small problems 
make fewer mistakes (and therefore learn faster) with the simpler 
languages of the small machines, provided that the problems to be 
solved are not too large for the small computers to handle conven- 
iently. If, however, the problems strain the memory size of the 
small computer to the point where the programmer must spend an 
appreciable part of his time figuring out a memory allocation 
scheme to fit the problems into the small computer, then the ad- 
vantage of language simplicity will be of little value. 
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The large computer will generally have superior diagnostics 
(tools for analyzing program errors) and more extensive subroutine 
facilities than the small computer. These aids are very much needed 
on the large, complex problems which big computers usually solve; 
large diagnostic and subroutine packages are not so essential on the 
smaller, more straight forward problems for which small com- 
puters are designed. 


Access and Turn Around Time 





In the area of access and turn around time, the small computer 
appears to come out ahead, since it will usually be located much 
closer to the department than would be a large computer which 
would be shared between departments, and the access of various 
problems to the small computer would be entirely controlled by the 
department, whereas on the iarge computer this access would be 
controlled by a priority system designed to accommodate the vary - 
ing needs of several different departments. In general, the farther 
removed the large computer is from the depariment which must use 
it, the greater will be the advantage in access and turn around time 
of the small machine. However, it should be borne in mind that this 
will only hold true for problems which are small enough to be con- 
veniently handled upon the small computer. 

Particularly during checkout of a new program, it is desirable 
for the programmer to be stationed near the computer so that he 
will be in a position to diagnose the results of his program and send 
it back to the computer as soon as possible. This usually results in 
sending the programmer to the computer for a sizeable period of 
time if the computer is remotely located from the programmer’s 
department, and adds to the cost of operating a remotely located 
computer installation. 

The complete control of access by one department means that 
the department, of course, will have unlimited access to the ma- 
chine, but that various individuals and groups within the department 
will still face a queuing or priority problem similar to that which 
the department would face in using a larger machine. The higher 
the load factor or the greater the percentage of time which the com- 
puter is actually in use out of a 24 hour day, the longer will be the 
turn around time and the more time will be lost by individuals wait- 
ing to make use of the computer. This will tend to be true on small 
systems as well as large systems, so that the access of any com- 
puter can only really be optimized if free time on the computer is 
available to take care of unexpected surges in demand. 
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Administration 


The administration problem as to whether to try to combine 
operations on a single large computer between various departments 
or to install separate smaller computers for the various depart- 
ments is usually the most important single factor in determining the 
feasibility of separate or joint operations. If departments fulfilling 
wholly different functions within a company are to share a single 
data processing installation, there must usually be some sort of a 
referee or administrative group set up to allocate computer time 
and decide the priority of various computer projects. If such a 
group does not exist, then the various departments will have to co- 
operate completely among themselves and if the departments tend to 
follow lines of authority on opposite sides of the company’s organi- 
zation chart (as for example in the case of Product Engineering and 
Accounting) then they may indeed have no well defined body to turn 
to to settle such priority and scheduling disputes as may be ex- 
pected to occur. For instance, the problem of what to do when a 
chief engineer insists on receiving test data by 3:00 p.m. and at the 
same time the hourly payroll is running three hours late has sel- 
dom been successfully resolved, except by the intervention of some- 
one in a position to arbitrate between the two departments involved. 
This sort of arbitration and compromise, of course, is not neces- 
sary if a small computer is used to serve each department. 

Since the big computer will require a certain amount of cooper - 
ation and indeed standardization between departments and since it 
has the capacity to handle very large problems, it also offers the 
potentiality of large scale data processing applications which will 
cut across departmental lines and offer efficiencies and economies 
which would not be possible in an independent solution of the prob- 
lems of any One department. This, of course, is one of the brightest 
goals of large scale data processing, to integrate the overall data 
processing requirements of entire companies. 

The price of such across-the-board integration of data process- 
ing is frequently a rather drastic reorganization of lines of depart- 
mental responsibility within the company and many companies have 
not yet found it practical to take a step of this type in an effort to 
integrate their data processing. When examined in this light, the 
choice of several small computers appears as a fairly conservative 
way to achieve reasonably limited goals, whereas the decision to 
pool the departments on a single large computer appears as a more 
radical way to achieve larger returns from the data processing in- 
vestment. 

The fundamental management policies of the company involved 
must necessarily play a large part in this decision. For instance, 
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if the company believes in a policy of complete decentralization and 
individual responsibility for its departments, it may tend to shy 
away from integrating the data processing of these departments into 
a single installation. Conversely, if the company is trying to cen- 
tralize management operations, it will probably find that integrated 
data processing on a single large machine will be a very powerful 
step in the direction of such centralized operations. 

In most small versus large computer decisions, administrative 
considerations outweigh the technical considerations of direct cost, 
throughput and access time and tend to establish a pattern of cen- 
tralized or decentralized data processing which reflects the man- 
agement philosophy of the company involved. 














wi 


A di 
p '(¢ 


rest 
of ¢ 
corm 
are 
terr 
proj 
rect 
cam 
ord 


is h 
ous 
reg 
tanc 
the 
roti 


is ¢ 
that 

















Parameters Needed for Stress Calculations 
on Cams 
By Robert Davies 


Research Laboratories Division 
General Motors Corporation 


Abstract. In this paper formulas are derived for the radius of 
curvature of a cam and for the distance along the cam and its fol- 
lower to the point of contact. If the lift p of a flat follower is the 
function p(@) of cam rotation 6, then 


radius of curvature = p”(6) + p(@), 


distance along follower = p’(@), 


0 
distance along cam = p’(0) + Jp (6) do. 
0 


A discussion is given of numerical means of calculating p’(@) and 
p”(@) when the function p(@) is defined by a table. 


Introduction. The stresses in a cam follower produced by its 
resting against the cam depend on the force of contact, on the radius 
of curvature of the cam, and on its width. The width and the static 
component of the force can be measured directly, but large errors 
are inherent in any scheme in which the radius of curvature is de- 
termined by constructing osculating circles on a drawing of the cam 
profile. It is more accurate to calculate the radius of curvature di- 
rectly from the cam-lift data or from the data used to grind the 
cam. Also, the acceleration of the cam follower must be known in 
order to calculate the dynamic force between the follower and cam. 

In examining the pitting of fatigued cams and cam followers, it 
is helpful to know what points on the surfaces are in contact at vari- 
ous cam rotations in order to determine the stress to which a pitted 
region is subjected. It is desirable to be able to calculate the dis- 
tance along the surface of the cam to the point of contact, as well as 
the distance across the face of the follower, as a function of cam 
rotation. 

Generally cams are prescribed by giving the amount of lift that 
is obtained as a function of cam rotation or by giving the amount 
that a cutter or grinding wheel is to be advanced toward the cam 
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Follower Face 
— x cos @ + y sin @ - p (6) = 0 


-_ 90 i 


/ P (6) \ 








Figure [The Equation for the Follower Face in Normal Form. 
axis as a function of cam rotation. The prescription is not equiva- 


lent to giving the rectangular or polar co-ordinates of the cam pro- 
file. 


equation of a straight line, that is, by 
x cos é@ + ysing -p(@) = 0, 


in which the x,y-axes are perpendicular to the axis of rotation and 
rotate with the cam and in which the origin is onthe axis of rotation. 
The perpendicular distance from the straight line to the origin is p, 
and @ is the angle between the perpendicular and the x-axis. 

If the cam follower is circular in cross section (a roller), or if 
the grinding wheel is circular, then the equation is 


[x - p(@) cos 6|° + [y - p(@) sin e}* = r’, 
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y Follower Face 
= cos 0, + y cin 6, ~p @) = 0 
90 
p (@) 
8, 
t 
v x 
8 
Cam 
x= p (8) cos ® - p' (8) sin 8 
y =p (@) sin 4 + p' (8) cos 8 
Figure 2. The Coordinates of a Point on the Face of the Cam 


and the Equation of the Follower. 


in which the same axes are used but in which Pp is the distance from 
the cam axis to the center of the roller, @ is the angle that a line 
between these points makes with the x-axis, and r is the radius of 
the roller. 

The cam profile is the envelope of all of these straight lines or 
circles. For a flat follower, the equation of the envelope is 


xcos 6 +ysin@ -p(@) =0. 
(1) -xsin 6 + ycosé@ - p’(6)= 0. 
The first equation represents the cam follower. The second is ob- 
tained from it by differentiating with respect to the parameter 64, the 
derivatives being indicated by primes. 

By eliminating @, the two equations can be reduced to give a sin- 
gle equation in x and y only (P is given as a function of @ and is 
therefore eliminated when @ is eliminated). The resulting equation 
is the equation of the cam profile. However, it is more convenient 
to work with Equation (1) in the parametric form: 
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x = p(@) cos 6 - p’(@) sin é, 


(2) y = p(@) sin 6 + p’(@) cos 6. 


In the case of a roller follower, the two equations are 


> 


[x - p(8) cos 6}* + [y - p(@)sine)*= r°, 


p() [x sin @ - y cos @] - p’(@) [x cos 6 + y sin 6} + p(9)p’(6) = O. 


Distance Along the Follower Face. The co-ordinates of the point 
of contact between the follower and the cam are the simultaneous 
solutions of the equations of the cam profile and of the follower. 
The equations of the cam profile are 





x = p(@) cos @ - p’(8) sin 6, 


y = p(@) sing + p’(@)cos @. 


in which @ is an unrestrained parameter. The equation of the cam 
follower is 


x cos §, + y sin#, - p(é,) = 0, 


in which 6, is the angle between the x-axis, a direction fixed with 
respect to the cam, and the normal to the cam follower face, a di- 
rection fixed in space. In other words, 6, is the angle of cam rota- 
tion. The simultaneous solutions of these equations are 


x = p(6,) cos 6, - p’(8,) sin 6,, 
(3) , 
y = p(6,) sin 6, + p’(6, ) cos 6, . 


It is evident that the point (x,y) on the cam profile, obtained from 
Equation 2 for a given value of the parameter 6, is the point in con- 
tact with the follower for the cam rotation @. 

The perpendicular from the axis of the cam to the follower face 
makes an angle 6, with the x-axis, and its length is (6, ). The co- 
ordinates of the intersection with the follower face are then 


x = p(6, ) cos 6, , 


(4) 
y = p(@,) sin 6, . 


This intersection is a convenient point from which to measure dis- 
tances along the follower face. The distance between the points 
given by Equations (3) and (4) is 
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which, therefore, is the distance along the follower face to the point 
of contact. At maximum lift, p’(@) is zero; consequently, distance is 
being measured from the point of contact at maximum lift. It is also 
the point closest to the cam axis. 


Radius of Curvature. In rectangular co-ordinates the radius of 
curvature R is 





2 2, 3/2 
. a7 oe 
x’y” a x"y’ 


for a curve in parametric form. The derivatives can be calculated 
directly from Equation (2). They are 


x’ = - [p(6) + p”(@)| sin 6, 

y’ = +[p(6) + p”(@)| cos 6; 
x” = -[p(@) + p”(@)| cos @ - [p’(@) + p”(6)) sine, 
y” = -[p(@) + p”(6)] sin 6 + [p’(@) + p”(@)| cose. 


Substitution in the equation for R gives 
R = p(@) + p” (9), 


which is the radius of curvature of the cam profile at (x,y), the point 
of contact for cam rotation 6. 


Distance Along Cam Surface. The arc length along a curve that 
is given in parametric form is 





© cients 
L= | yx’? + y’* do, 
1G, 
in which 6 is the value of the parameter at the point from which 


distance is measured. For the cam profile represented by Equa- 
tion (2), 


L = p’(@) - p’(@) + J p(é) dé 


is the distance between the points of contact for the positions @ and 
98, of the cam. Often, however, the distance along the cam is meas- 


0 


ured from the point of contact at maximum lift, and the rotation of 
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the cam is measured from the angle that gives maximum lift, so 
that 6, = 0 and p’(0) = 0. In this case, 


vg 


L = p’(@) + f p(e) dé 
0 


is the distance along the cam surface to the point of contact. 


Acceleration of the Cam Follower. The displacement p of the 
follower is 





p = p(9); 
the velocity and acceleration are therefore 


p = p’(a)8, 
p = p’(8)6 + p”(8)8’, 


in which dots denote derivatives with respect to time. 


For a cam that turns at uniform velocity, ¢ is zero and the ac- 
celeration of the cam follower is 


p = p”(9)6?. 


defined in such a way that it can be calculated nowhere except at 
limit points of the set over which the function is defined. Functions 
that are defined for every point in an interval may have a derivative 
in the interval; but functions defined over isolated sets of points — 
and these include all functions defined by tables — cannot have de- 
rivatives. Consequently, the function p(@), if it is defined by a table, 
must be extended to the intervals between the tabular entries before 
the derivatives p’(9) and p”(@) can be calculated. 

The cam consists of a sequence of minute flat facets when it is 
machined with a flat cutter. For such a cam it is possible to calcu- 
late the lift for all angles between the values in the table and then to 
calculate the corresponding derivatives. However, this method is 
not a very realistic one for interpolation because, in use, these 
facets soon loose their sharp edges and wear into a smoothly curved 
surface. 

Polynomials are frequently used to interpolate between tabular 
values when a smooth function is desired, since they are easy to 
evaluate and to differentiate. The interpolation polynomial for p(é@) 
can be found in the following way if the values of p(@) are given for 
equally spaced values of 6. Table 1 is a difference table; each entry 
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Table 1 - Example of Difference Table 


n p A A’ A A‘ A® A® 
J +23148 
+3801 
ff +26949 +396 
+4197 -14 
-1 +31146 +382 a 
+4579 -21 +1 
0 +35725 +361 -6 0 
+4940 -27 +1 
+1 +40665 +334 5 
+5274 -32 
+2 +45939 +302 
+5576 
+3 +51515 


is obtained from the column to the left of it by subtracting the entry 
immediately above it from the entry immediately below it. Table 2 
illustrates the notation that is generally used. For the differences 
in Table 2, 


Ap =A ‘p i 8 


m m+ m-1 
2 3 


The two left columns represent the value of n and p. The variable n 
is used in place of 6 in order to simplify some later equations. It 
is chosen to have the value zero in the middle of the table and to in- 
crease and decrease in increments of one. They are related by the 
linear transformation 


0 
(5) a*=—->5 @ = Rae + ¢, 
in which A@ is the increment in @ between entries in the table. 


There are no odd-order differences with subscript zero, but 
they can be introduced as the averages 
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Table 2 - General Difference Table 


n p rN A’ 4° a‘ a’ a’ 
on p_, 
4p _,3 
-2 p a p_ 
OP 43 Ap 
-1 p., ap, ' iii. 
Ap_. Ap, Ap , 
2 “3 - 
P 2 
0 Ps & p, Xp, A'p, 
Ap 1 ap, ap , 
2 "9 
+1 P., sp, Sp, 
4p. .3 Ap. 
+2 Pro ap,, i 
4p ot 
+3 Ps, 


The function 


n* 








> | 2 2.1 
p,(n) = Py + Bryn + ap, > + arp, ME) 4 wtp, BD) 
(6) 
5 n(n’-1)(n?-2? ) 2p n’ (n*-1)(n?-27) . .. (n? -(p-1)?) 
israeli. incite na (2p)! 


is obviously a polynomial in n and hence in @. Also, it takes on the 
2p + 1 tabulated values of p from n = -p to n = +p. The function is 
therefore a polynomial with which one can carry out the interpola- 
tion and differentiation. 

The derivatives are 








dp, : 5 3n-1  , «4n?-2n_ =, 55n*-15n?+4 

an 4° +Ap n+ Ap 31 +40, —ay +*&-—s7—* -- > 
dp, ’ . : 12n’-2 ' 20n* -30n 

ar OR TOR at On + 4aR Ct 
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dp} - I 1, 
dp, | 2 1 4 

(7) | : " @P «© Be «ess 
la} % ~ 12 °% 
| dP, A. 2 s 
ee m3 OO * OM 


The derivatives are very nearly equal to the corresponding tabular 
differences when the higher differences are small. For example, at 
n = 0 the first difference in Table 1 is 4759.5, and the relative error 
introduced by using it as the first derivative is TG = 0.0008. The 
relative error introduced by using 361 for the second derivative is 
aT = 0.001. This last formula, for calculating the derivative half 
way between the tabular entries, is given because of its greater 
convenience and more rapid convergence. It can, however, be de- 
rived in the same general fashion as the first equation. 

If the higher differences do not become small, Equation (7) can 
be used to calculate the correct derivative. However, tabular dif- 
ferences that eventually increase indicate that the table may contain 
errors. Table 3 shows how errors in a tabular value introduce 
larger errors in the differences, and how the errors become larger 
for the higher differences. The round-off errors in a table eventu- 
ally affect the differences sensibly, and the higher differences must 
be smoothed to remove these irregularities. 

Differentiation of Equation (5) gives 


dn 1 dn : 
dé ” Ag , dé? * ; 


The derivatives with respect to 6, which occur in the formulas for 
the radii of curvature and the distances, are 


(9) = # -292.i¢ 
f dé dn dé Aé dn’ 


‘ dp 1 d@ 
p(@) = SP = — SP. 
dé* Ae dn 


To calculate p’(@) and p”(@), the angle @ must be measured in ra- 
dians; therefore A? must be the interval in radians between tabular 
entries. 
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Table 3 - Difference Table with Error 


n ~p A A’ A’ a* A A® 
-3 p : 
Ap_ ot 
ot P_. . A’ p ; 
4p it A’p a 
-1 Pp, ‘ A’p_, +e ap 4 
Ap ; +€ A’p , -3¢ a’p , +106 
2 2 2 
0 pre Ap -2 r ap, +6€ Ap -20¢ 
0 0 
Ap. -€ Ap , +3« A’p. 1 -10¢ 
2 2 "3 4 2 
+1 p, ap, +e ap, -4e 
AP, 1 Ap. i-€ 
? 1s tls 
+2 _* ap, 
Ap 1 
+29 
+3 p., 


In order that Equation (7) may be used to calculate the deriva- 
tives, it must be assumed that the cam shape can be approximated 
by a polynomial. When this assumption is not valid, the derivatives 
may be in error. For a cam composed of flat facets, the lift can be 
calculated directly for all cam angles by a consideration of the geo- 
metrical configuration. The intermediate values so calculated are 
different from the values obtained from polynomial interpolation, 
and this polynomial interpolation is not valid for such a cam. It is 
not immediately obvious that the cam, as it wears into a smooth 
curve, comes any closer to a shape for which polynomial interpola- 
tion can be used. 


Calculation of fo(@) de. If the function p(@) is defined by a table, 
it must be extended to the intervals between tabular entries before 
the integral can be calculated. However, the value of the integral 
is less sensitive to the kind of interpolation used than are the de- 
rivatives, because integration is a smoothing process whereas dif- 
ferentiation is not. Generally polynomial interpolation is adequate, 
the only exception being the case in which the interpolated values of 
pP(8) ought to be systematically above (or below) a smooth curve 
through the tabular values. 
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The polynomial A, (n) given by Equation (6) is a suitable inter- 
polation function. Its integral is 
1 


2 2 3s 2 
a -. an 6 1 te Oe 
fo,@)dn=|pn+an > +ap, 5 +87 a 
“§ 
3n°_5n° " 1s: wT «2 
4 = —_— —_ 
OP. 513 % +24" - 5760 °% 


The integral is more conveniently calculated for values of n, or 6, 
that are half way between the values for which P(@) is tabulated. The 
increment between two successive values of the integral is very 
nearly the intervening value of (8), and one can calculate the inte- 
gral by adding the tabulated values of p(@), correcting the sum, if 
necessary, with the formula above. For the function shown in Ta- 
ble 1, the increase in the integral from n = -1/2 to n = +1/2 is 


361 . 17x6 
35725 + 94 * 5760 


= 35740. 
The term Je(e) dé is required in the formula for the distance 
along the cam surface. From Equation (5), 


dé = Aé dn 
and 
fo(@) dé = Aé JP, (n) dn. 


Once again, 6 must be in radians. The value of fe(e) dé for each 6 


is obtained by adding the corresponding Ag | % P, (n) dn to fe(e) de 
2 


for the previous @. 

The constant of integration is selected such that Jo(e) dé is zero 
for 6=0. It is necessary to evaluate p(@) dé at some point mid- 
way between the tabulated values of @ in order to start the iterative 
process for the calculation of subsequent values of Joe) dé. From 
Equation (6), 


€ Ad _€ 
iz) } = 4 = 
J Pl) ae Af JPp (n) dn = Ae 





+ A’P, 
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in which € is the fraction of the interval length A@ between zero and 
the next center ofan integral. If @ = 0 is a tabular value, then ¢ = 1/2 
and the integral is 


I 
3 406 


' - 1 1 2 
J ,Pt@) dé = Ad 5) p + 8 4p. ap. 


ee 

48 
At @ = 0, p(@) is a maximum and Ap is quite small. 

Summary. In calculating the stresses in a cam follower, one 


must know the radius of curvature of the cam and the acceleration 
of the follower. For a flat follower they are 


p 
radius of curvature = 7 + p(@), 


dé 
va l ti = fp dé . 
acceleration = ae dt 





in which p is the lift and @ is the angle of rotation of the cam in ra 
dians. 

When examining the areas of high stress in the cam and fol- 
lower, it is convenient to know the distances along their surfaces to 
the point of contact. For a flat follower they are 


distance along follower = se F 
dp ad 
distance along cam = ao * Jp(e) dé, 
0 


measured from the points in contact at maximum lift. 
In case P(@) is defined by a table, and polynomial interpolation is 
permissible, 


4 + A6 i — 7 < 
Je(e) dé = J (8) dé + Aé | dip + 94 4? - 5760 LU btp + +2]. 
dp _ 1 ee | 
dd ag | 4° - 3 AP 649 4° *---]> 
a Pee 
ao? 7 ast |S? - yp Se + | 
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The integral and first derivative are obtained for vaiues half way 

between the tabular entries, and the second derivative is evaluated 

at the entries (6 = 0 is an entry). The distance between entries is 
= Ad , 

Aé,in radians. The integral J p(@) dé is very nearly 9P(), and the 


summations are over the values after 6 = 0. 
























The Effect of Inclusions on the Fatigue 
Life of Steel 


By Leonard G. Johnson 
Research Laboratories Division 
General Motors Corporation 


INTRODUCTION 


It is generally admitted that non-metallic inclusions and other 
flaws in steel act as nuclei of stress concentration around which 
early fatigue failures can originate. This being so, it is important 
to be able to estimate what reduction in fatigue life can be expected 
for a given inclusicn content which has some representative stress 
concentration factor. By treating each volume element in a material 
as an independent component having a Weibull distribution of life, 
dependent upon its particular stress level, it is possible to make a 
theoretical prediction of the fatigue life for any specimen consisting 
of any number of such volume elements simply by multiplying to- 
gether the independent probabilities of failure of the separate ele- 
ments. In order to predict the effect of inclusions on fatigue life we 
need to know how many volume elements contain inclusions, and 
what the stresses are both in the elements containing inclusions and 
in those elements not containing inclusions. By means of these 
basic concepts a mathematical model of the effect of inclusions on 
fatigue iife can be constructed. Such a theoretical guide has been 
needed for a long time, and its construction will now be demon- 
strated. 


NOTATION 
v = Total stressed volume under consideration (grid of Fig. 1) = 1 
(consisting of k equally sized elements or components) 
k = Total number of components of volume under consideration 
V; = Volume of the i‘® component = : > : 
s = The nominal stress in the total volume under consideration 
s, = The stress in the it? component = A; 8 
A. = The stress concentration factor in the it? component 
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Fatigue Specimen 


c = The number of components of volume which contain inclusions 
(those with dots in Figure 1) 


v The sum of the volumes of all components having inclusions 
y 


f.% sy . fraction of the total stressed volume having inclusions 


POSTULATES OF THE MATHEMATICAL MODEL 





Postulate 1: The i‘? component (element of volume) has a life dis- 
tribution defined by the Weibull survivorship function [1]* 


where b is a positive constant having the same value for all com- 
ponents (known as the Weibull slope), 


Pp, (x) = exp | id (1) 


and where nj; = characteristic life of the i‘” component, i.e., the 
life for which the survival probability is e~* = .368. 


(pj (x) denotes the probability that the it® component will survive 
to a life of x.) 


Postulate 2: The relationship between the characteristic life of the 
i‘" component and the stress on the i‘" component is given by a for- 
mula of the type 


Kb 
~ P(s;) (2) 





ers in quare brackets refer to ref 
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where P(s;) is a monotone increasing function of s; (known as the 
potency of s;) (see [2]) 


(For pitting fatigue in bearnings, P(s;) is generally taken tobe a 
simple power function of s;, i.e., P(s;) = constant x s{" (m > 0)) 


Postulate 3: The components (elements of volume) comprising the 
specimen are statistically independent. 


DERIVATION OF THE LIFE PREDICTION EQUATION 





According to Postulate 1: 
; t 
p (x) = exp | - (x/n,)° | (3) 


Since according to Postulate 3 the components are statistically in- 
dependent we can arrive at the survival probability of the entire 
specimen by using the product rule |3| for independent survival 
probabilities. Thus, 


If p(x) = the probability that the entire assembly (specimen) sur- 
vives to life x 


It follows that 
p(x) = p(x) p (x) p (x) -****: P,, (x) (4) 


(k components) 
or p(x) = exp -x? (1 n? + 1/n? + 1/n? Ruin MO np) (5) 


Let n = the characteristic life of the assembly (specimen) 


Then p(x) = exp | -x>/a? | (6) 
", 1/i? = 1/n? + 1/2 + 1/n? + + 1/np (7) 
‘ : “ fig, Povevene /ny 
Let us denote by subscripts 1,2,3,....,c those c components having 
inclusions. 
Let n,,M,,M,,-..-,Mc be the respective characteristic lives of 


these c components. 


Let n = characteristic life of each component without inclusions 
(each at the same nominal stress s) 
(There are k-c such inclusion-free components) 


Then (7) becomes 


1/n> = 1/n> + 1/nP +....+ 1/n> + (k-c)/n> (8) 
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But, according to (2) 








b b 
i/n? = [P(s;)| and 1/n? - [P(s) | 
/ 1 k 1 = 
- qe = (P(s,)]° + (P(s, )}°+ siete + [P(s. > + (k-c)[P(s)]° 
' 6 kK 
or 1/q” = (PM 8)]" + [PA s)]°+.---- + [P(a_s)]” + (k-c)[P(s)]}° 
| k k 


Let \ = an “average” stress concentration factor such that 


(P(A, s)|> + [P(A,s) JP +...... + [P(A.s)]> = c¢ [P(As)]> 
Then 
-b e b b 
1/n’ = (c/k) [P(As)] + (1 - c/k) [P(s)] (9) 
or, 1/a® = £.[P(Xs)]” + (1 - f.) [P(s)]” (10) 


For the entire assembly (having unit stressed volume) we have (see 
Postulate 2) 


({n, = characteristic life of ) 

\ assembly at stress s 

and P(As) = 1/nj. (mds = characteristic life of 
assembly at stress A s) 


Thus, (10) reduces to 


-/aX, += (11) 


Equation (11) represents the fundamental life prediction equation 
for a material having the fraction f. of its total volume consisting of 
inclusions with a stress concentration factor 1. 


APPLICATION OF THE FUNDAMENTAL 
LIFE PREDICTION EQUATION 








Let us illustrate the use of the fundamental equation (11) by applying 
it to a bearing problem where the s--n relationship is 


n, =1/s™ (m > 0) (12) 
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nj, = 1/(A™ s™) (13) 
(11) becomes 


1/n” = f. Amb gmb 4, i f..) gmb = (1 - a 2 ,mb) gmb 


1 b 











n 
=b = i: = 4 a 
Fe 4 ole +0 2 - & 64, 
= ns 4 
we 8 4 wt. etek 
” 1 - 
n/n. (1 ¥ :. + £ AM) (14) 


Equation (14) gives us the ratio of the specimen life when the 
fraction f,. of the volume consists of inclusions of stress concentra- 


tion factor A to the life of an inclusion-free specimen. 


NUMERICAL EXAMPLE 





If the Weibull slope is b = 1, and if m = 12, it follows that if 1/1000 
of the volume of bearing material consists of inclusions of stress 
concentration factor A = 2 we can expect a fatigue life 


Ns 
n= 1 - .001 + .001 (212) = 19627 n . 





In other words, due to inclusions of stress concentration factor 2 in 
1/1000 of the volume, the life is reduced to less than 1/5 of what it 
would be for an inclusion-free bearing at the same nominal stress 
level. 


(The results for different stress concentration factors are shown 
graphically in Figure 2.) 


CONCLUSION 





By assuming Weibull life distributions of fixed slope and statisti- 
cally independent volume elements we have derived a simple life 
prediction equation which involves 


1. The fraction f. by volurme of material which contains inclu- 
sion. 


2. The stress concentration factor \ of the inclusions. 
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But, according to (2) 


1 


b b 
1/n? = {P(si)} and 1/n? - [P(s)} 








& 
‘ t/a? = [P(s,)]° + [P(s, )]°+ TTT + [P(s, yp + (k-c)[P(s) }° 
k k 
oe 1aP a (P(A, s)|° + (P(, s)]° Occden + (P(A.s)]° + (k-c)[P(s)]}° 
: k zs k 


Let \ = an “average” stress concentration factor such that 


(P(A, s)J> + [P(,s)JP+...... + [P(A.s)]® = ¢ [P(As)]> 
Then 
-b ; — b 
1/n = (c/k) [P(As)] + (1 - c/k) [P(s)] (9) 
or, 1/a” = £ [P(As)]” + (1 - f.) [P(s)]” (10) 


For the entire assembly (having unit stressed volume) we have (see 
Postulate 2) 


P(s) = 1/n, (n, = characteristic life et 
\ assembly at stress s 


and P(As) = 1/n\, (nj, = characteristic life of 
assembly at stress ts) 


Thus, (10) reduces to 


\= . « 
1/a® =f. /n®. ee (11) 

s 
Equation (11) represents the fundamental life prediction equation 
for a material having the fraction f. of its total volume consisting of 


inclusions with a stress concentration factor 1. 


APPLICATION OF THE FUNDAMENTAL 
LIFE PREDICTION EQUATION 








Let us illustrate the use of the fundamental equation (11) by applying 
it to a bearing problem where the s--n relationship is 


n, = 1/s™ (m > 0) (12) 


Ss 








ea aa oe 
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nj, = 1/(A™ s™) (13) 
(11) becomes 


P= 6 Pee +.) me Ok ot 











P nt 
Fe” 8 no he Oh Pm” 1 + 
_ ns 4 
~~ 6 O 6 «+t, +h 
: 1 2 
Wee" 4. f. + £,1™>)5 ated 


Equation (14) gives us the ratio of the specimen life when the 
fraction f, of the volume consists of inclusions of stress concentra- 
tion factor A to the life of an inclusion-free specimen. 


NUMERICAL EXAMPLE 





If the Weibull slope is b = 1, and if m = 12, it follows that if 1/1000 
of the volume of bearing material consists of inclusions of stress 
concentration factor A = 2 we can expect a fatigue life 


Ns 
1 - .001 + .001 (2'2) © 


n= 





19627 n . 


In other words, due to inclusions of stress concentration factor 2 in 
1/1000 of the volume, the life is reduced to less than 1/5 of what it 
would be for an inclusion-free bearing at the same nominal stress 
level. 


(The results for different stress concentration factors are shown 
graphically in Figure 2.) 


CONCLUSION 





By assuming Weibull life distributions of fixed slope and statisti- 
cally independent volume elements we have derived a simple life 
prediction equation which involves 


1. The fraction f. by volume of material which contains inclu- 
sion. 


2. The stress concentration factor \ of the inclusions. 





ion -Free 





1 


af 






| Bearing Fatigu 


onstant Stress 


Lili wti 


] 


l 


LEONARD G. 


Litiil 





JOHNSON 





| ecieedl 





001 

















LIFE 













INCLUSIONS ON FATIGUE OF STEEL 63 
Application and verification of this predicting equation requires that 
methods be devised whereby inclusion volumes and stress concen- 


tration factors due to inclusions can be accurately determined. 
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